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ABSTRACT 

Aims. We investigate statistical equilibrium of Cr in the atmospheres of late-type stars to show whether the systematic abundance discrepancy 
between Cr I and Cr II lines, as often encountered in the literature, is due to deviations from UTE. Furthermore, we attempt to interpret the 
NUTE trend of [Cr/Fe] with [Fe/H] using chemical evolution models for the solar neighborhood. 

Methods. NUTE calculations are performed for the model of Cr atom, comprising 340 levels and 6806 transitions in total. We make use of 
the quantum-mechanical photoionization cross-sections of Nahar (2009) and investigate sensitivity of the model to uncertain cross-sections for 
H I collisions. NUTE line formation is performed for the MAFAGS-ODF model atmospheres of the Sun and 10 metal-poor stars with -3.2 < 
[Fe/H] < -0.5, and abundances of Cr are derived by comparison of the synthetic and observed flux spectra. 

Results. We achieve good ionization equilibrium of Cr for the models with different stellar parameters, if inelastic collisions with H I atoms 
are neglected. The solar NUTE abundance based on Cr I lines is 5.74 dex with cr = 0.05 dex; it is ~ 0.1 higher than the UTE abundance. For the 
metal-poor stars, the NUTE abundance corrections to Cr I lines range from -1-0.3 to +0.5 dex. The resulting [Cr/Fe] ratio is roughly solar for the 
range of metallicities analyzed here, which is consistent with current views on production of these iron peak elements in supernovae. 
Conclusions. The tendency of Cr to become deficient with respect to Fe in metal-poor stars is an artifact due to neglect of NUTE effects in the 
line formation of Cr i, and it has no relation to peculiar physical conditions in the Galactic ISM or deficiencies of nucleosynthesis theory. 
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1. Introduction 

Abundance of chemical elements in the atmospheres of late- 
type stars is a key information in studies of Galactic chemical 
evolution (GCE). Combinations of different elements and vari- 
ation of their abundances with metallicity are commonly used 
to calibrate GCE models and to constrain poorly-known pa- 
rameters, like star formation history, initial mass function, and 
efficiency of mixing in the ISM. Since stellar yields are also a 
part of such models, the abundances can also probe the theories 
of stellar nucleosynthesis and evolution and highlight problems 
in them. 

The abundances are derived by means of high-resolution 
spectroscopy that implies modelling the observed stellar spec- 
trum. Whereas a poor quality of a spectrum introduces some 
random noise about a true value of the abundance, major sys- 
tematic errors result from an oversimplified treatment of radia- 
tion t ransfer and convection in stellar atmospheres (Asplund 
20051) . The assumption of local thermodynamic equilibrium 
(LTE), which is traditionally used to compute a spectrum in 



an attempt to avoid numerical difficulties with line formation, 
breaks down for the many species. Minority ions, which con- 
stitute at most few percent of the total element abundance, are 
particularly affected by NUTE conditions. As a consequence, 
various abundance indicators of the same element, like lines of 
different ionization stages or excitation potentials, often give 
discrepant results in ID UTE analysis. The main concern is that 
once observations compel us to restrict the analysis to a single 
indicator, e.g. due to a moderate signal-to-noise ratio or limited 
spectral range, spurious abundance trends with metallicity are 
unavoidable. 

For Cr, even-Z element of the Fe-group, only UTE cal- 
culations have been performed up to now. They revealed se- 
vere problems with modelling excitation and ionisation balance 
of Cr in the atmospheres of late-type stars. Systematic differ- 
ences of 0. 1 - 0.5 dex between abundances based on UTE fit- 
ting of the Cr i a nd Cr ii lines were reported for metal-poor gi- 
ants and dwarfs ([Johnson 2002 : Gratton et al. 2003 : Uaietal 



2008 : Bonifacio et a iJ l2009h . The discrepancies are smaller 
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in the atmospheres with larger metal content, a mounting to 
~ 0.1 de x for Galactic disk st a rs (p^ro chaska et al.'2000) and for 
the Sun dSobeck et af] 120071; lAsplund et al. 2009) . Following 
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Gratton & SnedenI (Il99lh . these offsets are usually attributed 
to the overionization of Cr i, a typical NLTE phenomenon af- 
fecting rninority atoms in stellar atmospheres. Alternatively, 
Bai et alj (120041) suggested uneven distribution of neutral and 
ionized Cr in different atmospheric layers, which depends on 
stellar parameters. 

As a result, two views on the evolution of Cr abundances 
in the Galaxy exist in the literature. The constant [Cr/Fe] 
with metallicity is derived from the LTE analysis of Cr ii lines 
JOratton & Sneden 1991; Bonifa cio et al ]!2009) that has a sim- 
ple interpretation in the theory of nucleosynthesis. Cr is formed 
with Fe in explosive Si-burning that occurs in supernova (SN) 
events, and the prodiiction ratio Cr/Fe is roughly solar in both 
SNe II and SNe la (Iciavtonll2003i) . Thus, using standard pre- 
scriptions for Cr nucleosynthesis, most of the G CE models re 



produce the fla t [Cr/Fe] trend with [Fe/H] (e .g. i Timmes et al 
1 9951; ISamlandlll 9981: iGoswami & Prantzoslfeom 



On the other hand, most LTE analysis of Cri lines 
JMc William et alJll995l:lc"avrel et al.l h004: Cohen et al. 200^ 
indicate that the [Cr/Fe] ratio steadily declines towards low- 
est metallicity. There is no simple explanation to this trend. 
Some chemical evolution studies tried to overcome the prob- 
lem by introducing correct i on factors to the theoretical SN 
yields (lFran9ois et al.l l2004t iHenrv et al. I l2010h . The assump- 
tion of erroneous supernova yields is not un realistic, given 
their sensitivity to the deta ils of explosion dlwamoto et al.l 
1999t iNakamura et al.l I1999 I). However, to reproduce an ap- 
proximate trend given by spectroscopic data arbitrary scaling 
factors to stellar yields are usually chosen, which lack any 
physical justification. Also, recent studies of metal-free mas- 
sive stars and their nucleosynthesis yields demonstrate that sub- 
solar Cr/Fe abundance ratios in very metal-poor stars can not 
be reproduced by any combinati on of SN II model parameters 
(lUmeda & Nomoto 2005; Heger & Woosl ev 2008), especially 
when other Fe-peak elements are taken into account. 

In this paper, we report NLTE abundances of Cr for the 
Sun and a sample of dwarfs and subgiants with -3.2 < [Fe/H] 
< -0.5. The atomic model for Cr and details about NLTE cal- 
culations are documented in Sect.|2] The statistical equilibrium 
of Cr under restriction of different stellar parameters in dis- 
cussed in Sect. |3] In Sect. |4] we present the NLTE analysis of 
solar Cr i and Cr 11 lines and derive the solar Cr abundance. In 
Sect. |5] we describe observed spectra, stellar parameters, and 
present Cr abundances for a sample of metal-poor stars. Finally, 
in Sect. |6] the revised abundance ratios [Cr/Fe] are compared 
with the trends predicted by GCE models and some implica- 
tions for the evolution of Cr in the Galaxy are discussed. 



2. The methods 

2.1. Statistical equilibrium calculations 

Restricted NLTE calculations for Cr are performed with the 
code DETAIL (Butler & GiddingsllT985l) . In the more recent 
version of the code, radiative transfer is based on the method 
of accelerated lambda iteration. 

The reference atomic model is constructed with 1 14 levels 
for Cr I and 225 levels for Cr 11, with energies of 0.01 eV and 1 . 1 
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Fig. 1. Grotrian diagram of the Cr i model atom. Lines represent 
allowed transitions included in the model atom. 



eV below the respective ionization limits, 6.77 eV and 16.49 
eV. The model is closed by the Cm: ground state. The typi- 
cal separation of fine structure components in Cr i and Cr 11 is 
less than 0.1 eV. Hence, we do not include fine structure in sta- 
tistical equilibrium (SE) calculations. Each term is represented 
by a single level with a total statistical weight of the term and 
energy weighted by statistical weights of fine structure compo- 
nents. Transitions are grouped. All transitions with oscillator 
strengths log gf < -8 and wavelengths A > 20000 A are ne- 
glected. As a result, the number of radiatively-allowed transi- 
tions is limited to 1590 and 5216 for Cri and Crn, respectively. 
Level excitation energies and oscillator strengths were taken 
from the Kurucz online databaseQ, supplemented by the exper- 
imental data of Murray (1992). The Grotrian diagram of the 
Cri model atom is shown in Fig. [1] The atomic data for some 
levels, including excitation energies Eexc, wavelengths A±i- and 
cross-sections cr^hi at ionization thresholds, are given in Table 

m 

The following atomic processes are taken into account 
in calculations of transition rates: radiative bound-bound and 
bound-free transitions, excitation and ionization by collisions 
with free electrons, and neutral hydrogen atoms. Scattering 
processes follow complete frequency redistribution. Recent 
calc ulations of ph otoionization from quintet and septet states of 
Cri (lNahadl2009h revealed resonances in cross-sections, which 
are due to the photoexcitation of the core leading to autoion- 
ization. The proper treatment of the latter process, e.g. tran- 
sitions to autoionizing states and ionization into excited core 
states, is not yet possible with the DETAIL code. Thus, we 
used partial cross-sections for ionization into the ground state 
of Cr II, which are lower than the total cross-sections at ener- 
gies above the first excited state of the core, A < 1500 A. This 
is a reasonable approximation for cool stellar atmospheres, be- 
cause temperatures are lower than those needed for core ex- 
citation. We have also shifted the cross-sections to the ob- 
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Fig. 2. Photoionization cross-section for the z^P° state of 
Cri: (black trace) quantum-mechanical calculations of Nahar 
(2009), (gray trace) hydrogenic approximation with effective 
principle quantum number «* = 1.88. 



served edge energies of Cr i states from the NlSld database, 
because the latter a re mor e accur ate than the calculated LS 
term energies from iNahar (l2009h . Photoionization from the 
other Cr i and Cr ii states is computed with a formula of Kramer 
( Menzel & Pek eris 1935) corrected for the ion charge (Rutten 
20031) and using the effective hydrogen-like principle quantum 
number n*. The use of effective quantum number increases the 
cross-sections at the ionization edge for low- lying levels and re- 
duces them for some levels of high excitation that is more real- 
istic than the original hydrogenic approximation. The quantum- 
mechanical cross-sections crQM(v) for the state z ^P° are com- 
pared with the hydrogenic approximation crhyd(v) in Fig. |2] 
The cross-section at the photoionization edge, as well as the 
background, are larger than crhyd(v). There are also strong reso- 
nances at energies, where the solar atmospheric UV flux is suf- 
ficiently strong to produce the overionization of the z ^P° state. 
Quantum-mechanical cross-sections for the other levels with 
excitation energies iiexc ~ 3 - 4 eV show a similar behaviour. 
We will show in Sect. [3] that these levels dominate ionization 
balance in Cri. 

At present, there are no experimental and accurate theo- 
retical data for ionization and excitation in Cri by electron 
impact, as well as by inelastic collisions with Hi. Thus, we 
rely on the commonly-used approximations. Cross-sections for 
continuum and allowed discrete transitions due to inelastic col- 
lisions with H I were computed with the formulae of D rawir 



19691) . This original recipe was developed for collisions be- 



tween equal hydrogen-like particles, hence we apply a scal- 
ing factor < S'h < 5 to the Drawin cross-sections. The fi- 
nal choice of Sn is discussed in Sect. 15.21 The rates of al- 
lowed and forbidden transitions due collisions with electrons 
are ca l culate d from the formulae of van Regemorter ( 1962) and 
Allenl(ll973h . respectively. The accuracy of the electron-impact 
cross-sections adopted in this work is not better than that for 
collisions with neutral atoms ( Mashonkina 1996). However, 
the excitation and ionization balance in Cri/Crii, as well as 



abundances of Cr determined from the lines of both ionization 
stages, depend only weakly on the exact treatment of inelastic 
collisions. Our calculations show that even in the solar at- 
mosphere ionization rates due to collisions with H i atoms dom- 
inate over the rates of ionization by electrons for the majority 
of the Cr levels, exceeding the latter by nearly three orders of 
magnitude for the uppermost levels. Moreover, collisional ion- 
ization rates dominate over collisional excitation rates. In the 
atmospheres of metal-poor stars, where the number density of 
free electrons is smaller than in the Sun, collisional excitation 
in Cri is fully controlled by neutral hydrogen. 

In addition to the reference atomic model of Cr, which 
is described above, we construct several test models. These 
models are used to show how the statistical equilibrium of Cr 
changes under variation of different atomic parameters, like 
size of the model atom, cross-sections for collisional and ra- 
diative transitions. 

2.2. Model atmospheres and spectrum synthesis 

To maintain consistency in the abundance calculations, we 
have decided to use the same type of model atmospheres, as 
employed for the derivation of stellar parameters (see Sect. 
|5]l, microturbulence velocities, and metallicites of the ob- 
jects investigated in this work. These are classical static ID 
plane-parallel models MAFAGS -PDF without chromospheres 
( Fuhrmann et al.|[T997 ). Line blanketing is treated with opacity 
distribution functions fromlKurucz ( 1992). Con vection is taken 
into account with the mixing-length theory (IBohm-Vi tense 



1958h and the mixing length is set to 0.5 pressure scale heights. 
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This value was chosen by Fuhrma nn et al.l (119931) to provide si- 
multaneously the best fitting of B ai mer line profiles in the solar 
flux spectrum with Tgff - 5780 K. iBarklem et al. (2002) also 
suggest that the best fit of Balmer line profiles can be achieved 
with a - 0.5. Stratifications of temperature and pressure in 
MAFAGS-ODF are similar to those given by other comparable 
models (see Fig. 15 in Grupp 2004). 

The abundances of Cr were computed by a method of 
spectrum synthesis with the code SIU, kindly provided by 
T. Gehren (private communication). Standard line broadening 
mechanisms, including radiation and quadratic Stark damping, 
were taken into account. The line half-widths due to elastic 
collisions with H i are calculate d using the cross-sections and 
velocity exponents tabulated by lAnstee & O'Maral (Il995b . In 
Table |2] they are given in terms of van der Waals damping con- 
stants log Ce, calculated for the temperature 6000 K. 

3. Statistical equilibrium of Cr 

The departures of atomic level populations n^^™ from their 
LTE values n^^^ are best understood from the inspection of de- 
parture coefficients, defined as /?, - n^^™/n^™. The departure 
coefficients for selected Cr i and Cr ii levels calculated for the 
solai- model atmosphere (reff = 5780 K, logg = 4.44, [Fe/H] = 
0, = 0.9 km s"') with different atomic models are shown as a 
function of continuum optical depth log Tc at 500 nm in Fig. [3] 
We will mainly discuss NLTE effects in Cr i, because the ma- 
jority of Crn levels remain very close to LTE for the range of 



4 



M. Bergemann & G. Cescutti: Cr in metal-poor stars 




log X (500 nm) log X (500 nm) 






-3 -2 -1 

log T (500 nm) 



1.2 r 



(h) 6000/4.2/-3.0 




-3 -2 -1 

log T (500 nm) 



Fig. 3. Departure coefficients of selected Cri and Cm levels as a function of optical depth at 500 nm. Panels (a) to (d) present 
results for the solar model atmosphere, (a): The reference model atom with a total of 374 levels constructed as described in Sect. 
12.11 Hydrogen collision rates are multiplied by 5h - 0.05. Photoionization cross-sections for Cri are taken from Nahar (2009). 
(b): Same as (a), but hydrogenic photoionization cross-sections are used for all levels (c): Same as (a), but 5h = 0. (d): The 
reduced model atom with a total of 108 levels including the Cm ground state, 5h = 0.05. (e) to (h): Same as (a), but the model 
atmospheres are taken from the grid with stellar parameters specified in each panel. Note that some levels on the plots are in 
thermal equilibrium with the other levels. 
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Table 1. Selected levels of Cri in the model atom. 
Photoionization cross-sections at the edge o-±r are adopted 
from Nahar (2009). 



Level 


Q 

5 


F 

^exc 


/111 


rr.i 
^ thr 


Level 


5 


F 

^exc 


/111 


^ thr 






eV 




Mb" 






eV 




Mb 


a'^S 


7 


0.00 


1832 


0.099 




5 


5.35 


8733 


0.316 


a^S 


5 


0.94 


2128 


0.316 


e^D 


25 


5.46 


9517 


3.504 


a^D 


25 


1.00 


2151 


0.327 


w^P" 


15 


5.48 


9645 


4.793 


a^G 


45 


2.54 


2936 


0.001 


v^P" 


15 


5.57 


10383 


5.26 


a^P 


15 


2.71 


3055 


2.133 




7 


5.66 


11195 


1.523 


z'P 


21 


2.90 


3209 


25.02 




5 


5.7 


11617 


2.613 




49 


3.15 


3431 


0.911 


/^D 


35 


5.8 


12768 


0.798 


z^F 


15 


3.32 


3599 


27.3 




15 


5.82 


13091 


29.07 




35 


3.42 


3708 


6.884 


v^F- 


35 


5.91 


14415 


20.47 




21 


3.45 


3741 


0.89 


g'D 


35 


5.92 


14554 


19.71 


V^P- 


15 


3.68 


4014 


3.756 


w^P° 


21 


5.92 


14565 


38.33 


Z^P 


35 


3.85 


4254 


0.255 


u^F 


35 


5.95 


15143 


0.028 


a^F 


35 


3.89 


4310 


0.547 


/^D 


25 


6.04 


16944 


4.298 


z'D" 


25 


4.17 


4779 


0.545 




7 


6.1 


18526 


2.968 


c'D 


25 


4.40 


5241 


19.14 




5 


6.12 


19033 


4.949 


e'S 


7 


4.57 


5656 


0.537 




10 


6.17 


20702 


10.83 


e'S 


5 


4.7 


5990 


0.903 




35 


6.22 


22667 


40.11 


x'P° 


15 


5.08 


7366 


16.38 


u^D° 


25 


6.28 


25238 


3.526 


y^F° 


35 


5.11 


7468 


0.091 


h^S 


5 


6.33 


28247 


7.905 


y^D° 


25 


5.15 


7687 


1.267 


t'D 


25 


6.45 


39357 


4.353 


z^H" 


55 


5.23 


8088 


0.009 


r^F 


35 


6.58 


68020 


0.129 


e'D 


35 


5.24 


8118 


6.231 


/I'D 


32 


6.62 


86092 


118.3 


x'P° 


21 


5.24 


8119 


27.37 


e''G 


63 


6.62 


86731 


151.0 




25 


5.29 


8378 


19.26 













" 1 Mb = 10-"*cm2 



Stellar parameters we are interested in. As seen in the Grotrian 
diagram (Fig. [U, Cr is a very complex atomic system because 
a partially filled 3d electron shell results in a highly complex 
term structure with many low-lying and closely-spac ed energy 



levels. Hence, in co ntrast to simpler atom s like Li dShi et al 



20071 Fig. 2) or Na (IBaumueller et al.lll998. Fig. 2), clear iso- 
lation of transitions driving departures from LTE is not trivial. 
We restrict the discussion of statistical equilibrium (SE) in Cr 
to few simple cases, which will be illustrated with a few levels 
characteristic for their depth dependency (Fig.|3]l. 

The NLTE effects in Cri are similar to that of the other 
minority atoms with complex configuration structure in cool 
stellar atmospheres, like Fe i, Mn i, and Co i. Deviations from 
LTE in Cr i develop in the layers where the mean intensity Jy 
exceeds the Planck function ByiJ^) over the bound-free edges 
of significantly-populated low-lying Cri levels. Analysis of ra- 
diative rates shows that the overionization is particularly strong 
from the odd levels with large photoionization cross-sections, 
e.g. z'^P" (^thr = 3209 A, Fig.©, z^F° (^thr = 3431 A), and 



z P° (/Ithr = 3599 A). Thus, their bj drop below unity akeady 
at logTc ~ -1-0.3 (Fig. [3^). Overionization affects other low- 
excitation levels with Egxc ~ 2 - 4 eV, although the net ion- 
ization rates are by orders of magnitude smaller then the rates 
from, e.g. z^P° level. The cross-section of the ground state a 
shows a number of autoionization resonances A < 1700 A, 



where solar fluxes are too low to produce significant ionization. 
The Cr i ground state is underpopulated via resonance transi- 
tion in the multiplet 4 (a - y ^P°), which is in detailed balance 
at logTc > -1.5, and via collisions with the lowest metastable 
states, which are overionized. A regular behaviour of depar- 
ture coefficients for the majority of Cr i levels at log > -2 
confirms the dominance of radiative overionization at these 
depths with some deviations from relative thermal equilibrium 
between low levels (compare levels z ^F° and a in Fig. [3^) re- 
flecting the non-hydrogenic character of their photoionization 
cross-sections. We have also performed SE calculations with 
photoionization cross-sections for all Cri levels derived from 
hydrogenic approximation. In this case, departure coefficients 
are homogeneously distributed up to the depths logTc ~ -1.5 
(Fig.[3j3); bi is a nearly monotonic function of the level ioniza- 
tion potential. 

Line transitions influence Cri excitation balance in the 
outer layers, log < -2, where departure coefficients of the 
levels markedly deviate from unity and from each other (Fig. 
[3^). Pumping by superthermal radiation with Jy > ByiT^) be- 
comes important in the layers, where optical depth in the wings 
of strong low-excitation lines drops below unity. Detailed bal- 
ance in the transitions does not hold and their upper levels are 
overpopulated. As soon as line cores form, Tq < 1, spontaneous 
transitions depopulate the upper levels and their bi steeply de- 
crease. The spectrum of Cr i is represented by a large number of 
lines in the near-UV, which are subject to these non-equilibrium 
excitation effects. The transitions in the multiplet 4, which sus- 
tain thermal equilibrium between a ^S and y ^P° at log > -1, 
now go out of the detailed balance and there is a small net radia- 
tive absorption leading to the overpopulation of y^P" (Fig. [3^). 

At log Tc 3, photon pumping ceases, because in the cores 

of the resonance lines rj, < 1, hence the departure coeficient 
of the y''¥° level drops. The levels with Eexc ~ 2.9 - 3.4 eV 
are radiatively connected with the uppermost levels, which are 
separated by an energy gap < 0.5 eV from the continuum, e.g. 
transitions y ^P° - e ^F with A ~ 3880 A. Net radiative rates in 
these transitions are positive at log > -2, but overpopulation 
of the upper levels does not occur because they are strongly 
coupled to the fully thermalized Cm ground state. This cou- 
pling is maintained by collisions with Hi, which exceed the 
respective rates of ionization due to collisions with electrons 
by few orders of magnitude. There is also a sequence of spon- 
taneous de-excitations through the densely packed upper levels 
with Eexc ~ 5 - 6 eV, which is driven by Jy < By(Tg) at the 
corresponding frequencies. 

Fig. [3]: demonstrates the behaviour of Cri level popula- 
tions, when inelastic collisions with H i are neglected, Sh = 0. 
Underpopulation of the levels at all optical depths is amplified, 
also the uppermost levels, like e ^F, decouple from the contin- 
uum already at logTc ~ 0. Large differences with the refer- 
ence model (Fig. [3^), computed with Sh = 0.05, are seen at 
logTc < -L5, where electronic collisions are also not effec- 
tive. Hence, we expect significant changes in the opacity of the 
cores of strong Cr i lines under NLTE. Also at log Tc > -2, de- 
parture coefficients of high-excitation levels are different from 
the case of Sh = 0.05, thus affecting the source functions of 
intermediate-strength lines. 
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Similarly, removal of the high-excitation Cri levels leads 
to slightly increased deviations from LTE in Cri, although the 
variation of departure coefficients with depth at -2 < log Tc < 
is not different from the reference complete atomic model. 
Results for the reduced model atom, constructed with 108 lev- 
els of Cri with excitation energy of the highest level 6.76 eV 
and closed by the Cr ii state, are shown in Fig. |3}l- This model 
is devoid of some doubly-excited Cm states below the 1-st 
ionization threshold and only transitions with \og gf > -1 
are included. Note the amplified underpopulation of the low 
and intermediate-excitation levels compared to the reference 
model. This result can be easily understood as due to the re- 
duced radiative and collisional interaction of the Cr i levels with 
each other and with the continuum. Less electrons recombine 
and de-excite to the lower levels via weak infrared lines, where 
Jv < B,(T,). 

Statistical equilibrium of Cr in the atmospheres of cool 
subdwarfs and subgiants is established by radiative processes. 
NLTE effects on the levels of Cr i and Cr ii are amplified com- 
pared to the solar case. The main reason is low abundances of 
metals, which supply free electrons and produce line blanket- 
ing in the short-wave part of a spectrum. Hence, increased UV 
fluxes at bound-free edges of low-excitation Cr i levels lead to 
their strong underpopulation, and collisional coupling between 
the levels is very weak due to deficient electrons. As an exam- 
ple, we consider a metal -poor turnoff star with [Fe/H] - -2.4 
(Fig. [3^). Ionization balance is dominated by radiative transi- 
tions from the low-excitation Cri levels. At -1 < log < 0.2, 
these levels are depopulated by overionization and by optical 
pumping in strong resonance lines, connecting the levels of 
a^S andy^P" terms. At logTc < -1, photon pumping ceases 
because the lines become optically thin, and spontaneous de- 
excitations maintain relatively constant /?, of the a ^S state in 
the outer layers. An interesting result is the presence of small 
deviations from LTE for intermediate-excitation Cr ii levels. As 
seen on the Fig. [3^, the odd level z'^D° with excitation en- 
ergy 6.2 eV is slightly overpopulated by optical pumping at 
-2 < log Tc < 0. But at smaller depths, photon losses in transi- 
tions of the multiplet b*F - z'^D° result in a depopulation of the 
upper level. We use the lines of this multiplet in the abundance 
analysis. 

Effective temperature and gravity do not affect distribution 
of atomic level populations at solar metallicity, but they be- 
come important with decreasing [Fe/H]. The results for the 
model with T^ff = 6400 K, logg = 4.2, and [Fe/H]= (Fig. 
[3]0 are almost indistinguishable from those obtained with the 
solar model atmosphere at optical depths logTc > -3. This is 
stipulated by increased collisional interaction among the Cri 
levels with fiexc > 3 eV. Although the model flux maximum is 
shifted to shorter wavelengths, ground state ionization is still 
not efficient because the edge cross-section of the a ^S state 
is very low. In the cool model atmospheres, deviations from 
LTE depend on the stellar gravity, and the effect is most pro- 
nounced at low metallicity. At [Fe/H] = -3, departures from 
LTE for Cr i levels are stronger in the model with Teff - 5000 
K and logg = 2.6 then in the model with Teff = 6000 K 
and logg - 4.2 (Fig. |3^, h). This may account for a system- 
atic difference between metal-poor giants and dwarfs found by 



Lai et al.l (120081) and iBonifacio et al.l (120091) (see discussion in 
Sect.ESl. 

The atomic model of Cr used in this work is not compu- 
tationally tractable in full NLTE calculations with 3D model 
atmospheres (see e.g. lBotnen&Carlssonlll999l for Ca), which 
require reduction of the model. Our analysis suggests that it is 
possible to construct a simpler model of Cr atom, which in- 
herits the basic properties of the complete model and has sim- 
ilar performance under restriction of different stellar parame- 
ters. In particular, removal of high-excitation levels and nu- 
merous weak transitions in the Cri atom, whose main effect 
is to provide stronger coupling of levels, can be compensated 
by increasing collisional interaction between them. The care- 
fully chosen scaling factor to inelastic H i collisions makes up 
for the missing transitions: for Cr i, the reduced model with 108 
levels and Sh - 0.15 gives a similar description of the statisti- 
cal equilibrium in Cr for solar-type stars to the complete refer- 
ence model with 340 levels and - 0.05. However, there are 
two important concerns. First, the fact that the simple model 
atom performs well in ID does not guarantee that accurate re- 
sults are obtained with the same atomic m odel with 3D con- 
vective model atmospheres ( Asplundl 2005h . Also, substituting 
the multitude of levels and transitions in Cr i by increased ef- 
ficiency of collisions with H i (although both seem to produce 
the same effect on level populations) has certainly no physical 
justification. In fact, our result even suggests that large scaling 
factors to the Drawin's formula, as sometimes encountered in 
the literature, may stem from the missing atomic data in SE 
calculations. 

4. Solar lines of Cr i and Cr ii 

All lines selected for the solar abundance analysis are given in 
Tabled We used the MAFAGS-ODF model atmosphere (Sect. 
IS with solar parameters. The profiles are broadened by the 
solar rotation Vi-ot.o = 1.8 km s"' , microturbulence velocity 
= 0.9 km s"' , and by a radial-tangential macroturbulence 
velocity Vmac - 2.5 . . .4 km s"'. Radial and tangential compo- 
nents of Vmac are assumed to be equal, and velocity distribu- 
tion for each component is Gaussian. V^ac is allowed to vary 
with the line strength and depth of formation. Th e comparison 



spectrum was taken from the Solar Flux Atlas of iKurucz et al 
(1984). 

Oscillator strengths for the Cri lines are taken from two 
sources. Most of the data are from ISobeck et al.l (l2007h . who 
determined branching ratios of transitions by Fourier trans- 
form spectroscopy. The accuracy of logg/ values is within 
5 - 10% for majority of transitions. In the absence of data 
from Sobeck e t al., w e used transition probabilities from 
Blackwell et al. (1986) measured with the absorption tech- 
nique. iFuhr et al.. (J988.) ascribe 10% accuracy to these val- 
ues. Oscillator strengths fo r the majority of Cm lines are taken 
from Nils son et al.l (12006 ). Branching ratios were measured 



with Fourier transform spectrometer and combined with life- 
times from laser-induced fluorescence experiment. The uncer- 
tainty of absolute oscillator strengths is 12 - 16% for the tran- 
sitions used in this work, except for the Cm line at 4592 A 
with an uncertainty of 37%. For the lines at 4555 and 4812 A, 
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Table 2. Lines of Cri and Cm selected for solar and stellar 
abundance calculations. Wavelengths A and excitation energies 
of the lower levels of the transitions E\o^ are taken from NIST 
database. The multiplet is specified in the 3-d column. Solar 
logigfe) values for the lines with an asterisk in the wavelength 
entry can not be reliably computed due to severe blending. 
These lines are not used in determination of the solar Cr abun- 
dance. 
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we used the logg/'s from IWuiec & Wenigen (Il98ll) . and the 
Cr II lines of multiplet 43 were calcul ated with log.g/'s derived 
from the solar spectrum analysis by Kostvk & Orlova 
The accurac y of log^/ value s from Wuiec & Wenig er (1981) 



and jKostvk & Orlovai d 19831) was later revised by iFuhr et al 
dlMSbto be of the order of 50% 



The LTE abundance of Cr determined from the Cri lines 
is logs = 5.66 + 0.04 dex, where the uncertainty is one stan- 
dard deviation. Our LTE abundance is consistent with the me- 
teoritic value, 5.63 + 0.01 dejJl Sobeck et al. (2007) calculated 
loge = 5.64 + 0.04 dex from equivalent widths of Cri lines in 
the disk-center intensity spectrum. The agreement of our LTE 
abundance and that of Sobeck et al. is surprising, since these 
authors used a semi-empirical Holweger-Miiller model atmo- 
sphere and assumed a lower value of the microturbulence ve- 
locity, = 0.8 km s ' . The abundances of Cr based on the LTE 
analysis of Cm lines are similar in our study and in Sobeck et 
al., 5.81 + 0.13 dex and 5.77 dex with cr = 0.13 dex, respec- 
tively. 

The LTE analysis of Cr lines revealed two problems. First, 
there is an abundance discrepancy of ~ 0.15 dex between Cri 
and Cm lines, which also appears in the study of Sobeck et al. 
The line-to-Une abundance scatter for Cr 11 is also large, cr - 
0.13 dex. Second, there is a large spread of abundances within 
the multiplet 18: profile fitting of the Cri lines with equivalent 
widths Wa > 60 mA requires systematically higher abundances 
by ~ 0.1 dex compared to weaker lines. A similar discrepancy 
was noted by Blackwell et al. (1987), who explained this with 
non-thermal excitation in the multiplet 18. Sobeck et al. did not 
support the Blackwell's conclusion, although a 0.1 dex line- 
to-line scatter in this multiplet is also present in their results. 
We show in the next paragraph, whether either of these two 
problems can be solved with NLTE. 

The NLTE line formation is based on the departure coeffi- 
cients bi computed with the reference model of Cr atom (Sect. 
|2]i. We performed test calculations for two scaling factors to 
inelastic collisions with hydrogen, S'h = and 0.05. In both 
cases, NLTE corrections to abundances derived from Cri lines 
are positive and range from 0.05 to 0. 1 dex, depending on the 
line strength and the excitation potential of the lower level. 
All Cr I lines computed under NLTE with Sh = fit the ob- 
served spectrum with the Cr abundance log s - 5.74+0.05 dex, 
whereas a scaling factor Sh - 0.05 leads to log s = 5.7 + 0.04 
dex. The NLTE abundance corrections to the Cm lines are 
small and negative. The results for 5h = and Sh - 0.05 
are equal, logs - 5.79 + 0.12 dex. In contrast to LTE, the 
difference between both ionization stages is now fairly small, 
log e(Cr II) - log e(Cr i) = 0.05 dex for Sn - 0, and it is within 
the combined errors of both values. At present, we can not de- 
termine what causes the discrepancy between the abundance 
of Cr in meteorites and our NLTE ID abundance for the solar 
photosphere. The accuracy claimed for oscillator strength s of 
Cr II and Cr i t ransiti ons measured bv lNilsson et al. ( 2006h and 



Sobeck et alj (120071) . respectively, is very high. There is evi- 



^ The Cr abundance in Cl-chondrites from iLodders et al.l ll2009l) 
was renormalized to the photospheric Si abundance of Ishiet^ 
(2008) . log £si.o = 7.52 dex 
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dence that lines of a similar atom Fe i are affected by convective 
surface inhomogeneities, also NLTE effects on Fe i lines are 
different in ID and in 3D radiative transfer calculations (see 
below). For Cr, such detailed investigations are not yet avail- 
able. 

Our reference NLTE model satisfies ionization equilibrium 
of Cr, but it does not solve the problem of overestimated abun- 
dances from several lines of multiplet 18. The solar NLTE 
abundance correction^ Anlte for all lines of this multiplet are 
equal and they scale with the Su parameter. The behaviour of 
line profiles under NLTE can be understood from the analysis 
of level departure coefficients bj at the depths of line formation. 
The formation of weak lines with Wa < 60 mA is confined to 
- 1 < log Tc < 0, where bj = bj, so line source functions are 
thermal, = By. Hence, line intensity profiles reflect the pro- 
file of an absorption coefficient, Ky ~ b,. Since fe, < 1, the lines 
are weaker relative to their LTE strengths and Anlte is posi- 
tive. Radiation in stronger lines of the multiplet 18 comes from 
the depths -4 < log < 0. The line cores have lower inten- 
sities under NLTE due to depleted source functions, ^ij < By, 
at logTc < -2. But their wings, which dominate the total line 
strength, are formed at the depths where fe, < 1 due to overi- 
onization. The net effect on the line profile is that the NLTE 
abundance correction is positive. For all lines of multiplet 18, 
we find Anlte ~ +0.1 dex and -i-0.05 dex for 5h = and 
Sn - 0.05, respectively. The result is incompatible with the 
suggestion of lBlackwell et al. ( 1987 ). since the latter would re- 
quire NLTE effects of different magnitude for all lines of multi- 
plet 18. Moreover, in our 'extreme' NLTE model with = 0, 
the net NLTE effect on the profiles of stronger lines must be 
zero to eliminate the abundance scatter If this problem is due to 
a deficiency of the atomic model, our main concern is the crude 
approximation used for inelastic collisions with H l Compared 
to the Drawin's formulae, ab initio quantum-mechanical cal- 
culations dBelvaev et alJll999HBelvaev & Barklemll2003h pre- 
dict significantly lower collision rates for certain transitions of 
simple alkali atoms, and they show that, in addition to exci- 
tation, other effects like ion-pair formation become important. 
Thus, in principle, any rescaling of the Drawin's cross-sections, 
e.g. using variable 5h for levels of different exci tation energies 
as exemplified by iBaumueller & Gehren (Il99d) . does not im- 
prove the reliability of results. Variation of Sn does not reduce 
the scatter among the Cr n lines, because the NLTE effects are 
very small even for Sa - 0. For the same reason, employment 
of an LTE versus a NLTE approach does not decrease the scat- 
ter. 

The abundance anomaly among some Cri and Cm lines 
may reflect shortcomings of the ID static mixing-length model 
atmospheres. Discrepant Cri and Cm lines with 70 < < 
120 mA are very sensitive to a variation of the microturbulence 
parameter that was also demonstrated for Fe i lines with simi- 
lar equivalent widths (Gehren et al. 2001b). For A^t - ±0.2 km 
s ' with respect to the reference value 0.9 km s the abun- 
dances derived from the Cri lines of multiplet 18 change by 




5247.5 



5247.6 
Wavelength (A) 



5247.7 



The difference in abundances required to fit LTE and NLTE 

profiles is referred to as the NLTE abundance correction, Anlte = 
loge^^-logfiLTE 



Fig. 4. NLTE profiles of the Cr i fine at 5247 A. Different com- 
binations of microturbulent velocities and resulting Cr abun- 
dances log E are shown. The line shows a convective asymme- 
try. 



roughly +0. 1 dex, and any combination of and log e leads to 
an equally good profile fit (Fig.lH). Using the depth- dependent 
microturbulence as suggested bv iHolwegei ( 1967 ) does not 
eliminate the abundance scatter among Cri lines, although the 
Cr II lines are strengthened requiring 0. 1 - 0. 15 dex lower abun- 
dances. 

Clearly, our approach of solving line formation in ID ho- 
mogeneous models with a two-component gaussian velocity 
field is simple com pared to multidimen s ional radiative trans- 
fer cal culations (e.g. lAsplund et al ]l2000h .l ISteffen & Holwegen 
(l2002h lave shown that abundance corrections due to photo- 
spheric temperature fluctuations vary in sign and magnitude 
for lines of different E\o^, A, and Wa- For saturated lines, some 
effect due to hydrodynamic velocity fields are also predicted. 
Besides, there is evidence that NLTE effects on atomic level 
populations are ampli fied in the presence of convective sur face 
inhomogeneities (e.g. Shchukina & Truiillo Buenol 2001 . for 
the Sun). Neutral iron is strongly overionized above hot gran- 
ules compared to cooler intergranular regions with the max- 
imum effect on the opacity of low-excitation Fe i lines with 
Eiow < 2 eV. Deviations of the source function from By{Ts) 
are important for higher-excitation lines in the intergranular re- 
gions. Shchukina & Trujillo Bueno also find that strong lines 
of any excitation potential formed above log 2.5 are sen- 
sitive to the variation of collision rates. Comparable effects can 
be expected for the lines of Cr i, hence a s ystematic study of 3D 
effects in Cr is necessary to confirm this. lAsplund et al. [ (l2009h 
computed LTE line formation for Cr i and Cr ii with a hydrody- 
namical model atmosphere. Applying NLTE abundance correc- 
tions derived from our ID modelling, Asplund et al. obtained 
the solar Cr abundance 5.64 ± 0.04 dex, which is almost equal 
to our LTE ID solar abundance determined from the Cri lines. 
Another study, which provides an order-of-magnitude estimate 
of 3D effects on the Cr line f' o rma tion in very metal-poor at- 
mospheres, is iBonifacio et al. ( l2009i) . We briefly discuss this 
paper in Sect 
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Finally, we note that all problematic Cr ii lines are located 
in the spectral region 4500 - 4800 A, where the continuum 
level is uncertain due to severe blending. Hence, the abundance 
discrepancy may be, at least in part, removed by adjusting the 
local continuum. On the other side, this procedure could also 
affect Cr i lines located in this spectral region. 

5. Abundances of Cr in metal-poor stars 

5.1. Observations and stellar parameters 

The list of stars along with their parameters is given in Table 
[3] Most of the stars were observed by T. Gehren and collabora- 
tors with UVES spectrograph at the VLT UT2 on the Paranal, 
Chile, in 2001, and/or with the FOCES echelle spectrograph 
mounted at the 2.2m telescope of the CAHA observatory on 
Calar Alto. The observations of HD 84937 were taken from the 



UVESPOP survey dBagnulo et al..2003.) . The UVES spectrum 
for G 64-12 was taken from ESO/ST-ECF Science Archive 
Facility (FID 67.D-0554(A)). With the spectral resolution of 
A/AA ~ 40000 - 60000, it is not possible to distinguish the 
contribution of Vm^c and Viot to a line profile. Hence, the pro- 
files are convolved with a Gaussian with the full width at half 
maximum ~ 2.8 - 4 km s"'. 

The sample includes nine objects from our recent studies 



jBergemann & GehrenI l2008t iBergemann et al. I l2009l) . where 



the details on observational material and derivation of stellar 
parameters can be found. In addition, we include here one thin 
disk star, HD 123710. Parameters of all these stars are deter- 



mined bv lGehren et all (i2004., 2006) using the same MAFAGS- 
ODF model atmospheres, as employed in this work. The ef- 
fective temperatures were obtained from fitting the Balmer 
line profiles under LTE. The spectroscopic temperatures are in 
agreem ent with the Tjff's determined by the method of Inf rared 
Fluxes / Alonso et al.]|l 99611 1 9991 ICasagrande et al.ll2010l) . The 
off'set with Teff's from the last reference is ~ 20 K with an 
nns error 23 K. Surface gravities are based on Hipparcos par- 
allaxes. 

Metallicities and microturbulence parameters were deter- 
mined from LTE fitting of Fe ii lines by requiring a zero slope 
of Fe abundances versus line equivalent widths. The validity of 
LTE approach in ID line formation calculations for Feii was 
demonstrated in very accurate analyses o f statistical equilib- 
rium ofTeJntheSmM Gg^ 200 la b) and in metal-poor 
stars (iMashonkina et al ■I l2010h . These studies employ the same 
typ e of atmospheric models, a s we use in the current work. 

Gehren et all (l2004 EoO ^ estimate the errors to be 100 K 



for Teff, 0.05 dex for \ogg, 0.05 dex for [Fe/H], and 0.2 km/s 
for For the halo and thick disk stars, the model atmospheres 
were computed with a total (T-element enhancement of 0.4 dex. 

5.2. [Cr/Fe] abundance ratios 

The abundances of Cr in metal-poor stars are calculated strictly 
relative to the Sun. Any abundance estimate derived from a 
single line in a spectrum of a metal-poor star is referred to 
that from the corresponding solar line, which excludes the use 
of absolute oscillator strengths and of the average solar abun- 



dance. We derive a diff'erential element abundance for each de- 
tected Cr line in a spectrum of a star according to 

[El/H]=logfe/£)*-log(g/ef 

and then we average over all these lines for a given ionization 
stage. Most of the Cr lines analyzed in the solar spectrum be- 
come very weak at low metallicities, so typically we use from 
3 to 10 lines for each ionization stage. 

Fig. |5] displays the diff'erence between [Cr/Fe] ratios com- 
puted from the lines of two ionization stages, [Cru/Cri] - 
[Crn/Fe] - [Cri/Fe], as a function of stellar effective tempera- 
ture and metallicity. NLTE and LTE ratios are shown with filled 
and open symbols, respectively. NLTE abundances are based 
on SE calculations with the atomic model of Cr with Sh - 0. 
There is an offset of ~ 0.2 - 0.3 dex between Cri and Crii 
lines under LTE, which increases with decreasing [Fe/H] (Fig. 
|5] bottom). This reflects the fact that the main stellar parameter 
that controls deviations from LTE in Cr is the metal abundance. 
Since the NLTE effects of gravity and effective temperature are 
smaller than that of [Fe/H], we do not see a clear trend with 
Teff in Fig.|5](top), where no pre-selection according to log g or 
[Fe/H] was made. A discrepan c y of a similar ma gnitude and 



sign was repor t ed by Ijohnsoni (l2002h . iLai et al.l (l2008l) . and 



Bonifacio et al. ( 2009h . All studies agree that for giants with 



Tes < 5000 the difference between two ionization stages is 
~ 0.4 dex. For hotter and higher-gravity models, characteristic 
of dwarfs, Lai et al. and Bonifacio et al. derive smaller offsets, 
[Cru/Cri] ~ 0.2 dex, in good agreement with our LTE calcu- 
lations. These findings are consistent with differential NLTE 
effects in Cr, which at low metallicity are larger for low log g 
and low Teff models corresponding to giants (Sect. [3]). 

Our main result is that the offset between Cr i and Cr ii for 
any combination of stellar parameters investigated here disap- 
pears when we use atomic level populations computed with the 
reference atomic model with 5h = (Fig. |5]). This model al- 
lows for extreme NLTE effects on Cr i due to neglect of ther- 
malizing Hi collisions. For Sh = 0.05, the abundances based 
on the Cr i lines decrease by < 0.05 dex. Although these results 
are close to the abundances computed with Sh = 0, the agree- 
ment between two ionization stages is not so good, [Cru/Cri] 
~ 0.04 . . . 0.07 dex, and the offset becomes even larger with 
increasing Sh- 

The average Cr abundances for each star and their stan- 
dard deviations are given in Table [3] Table 4 with individ- 
ual abundances for each line of Cri and Cm is available 
in the electronic edition of A&A. [Cri/Fe] and [Cru/Fe] re- 
fer to the analysis of Cri and Cm lines, respectively. NLTE 
abundances given in this table were computed with Sh - 0. 
In the spectra of BD-4°3208 and G 64-12 no lines of Cm 
were detected, and [Cru/Fe] value for HD 140283 is based on 
one Cm line. In Fig. |6] [Cr/Fe] ratios computed under NLTE 
and LTE are plotted as a function of the stellar iron abun- 
dance, with [Fe/H] based on the LTE analysis of Feii lines. 
The mean NLTE [Cr/Fe] ratio in stars with subsolar metal- 
licity is ([Cr/Fe]) = with the dispersion cr - 0.06 dex. 
Assuming LTE, we derive ([Cr/Fe]) = -0.21 with cr = 0.11 
dex. The mean LTE value ([Cr/Fe]) and the slope of [Cr/Fe] 
with [Fe/H] agree with other comparison studies, which refer 
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Table 3. Stellar parameters for the selected sample. NLTE and LTE abundances of Cr computed using the lines of Cri and Crii 
are given as [Cri/Fe] and [Crn/Fe], respectively. The standard deviations of these values are quoted as errors. The abundances of 
Fe are based on Fe ii lines. Nci i and Ncr ii are the number of Cr i and Cr ii lines used for each star. See text for further discussion 
of the data. 
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[Mg/Fe] 
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HD 19445 
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1.5 


-1.96 


7 


3 


-0.07 = 


t0.05 


-0.3 ± 0.06 


-0.09 + 0.01 
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1.7 
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10 


2 
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0.24 
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14 


4 
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-0.04 ± 0.03 
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-0.02 + 0.02 


0.24 


HD 148816 
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4.07 


1.2 


-0.78 


21 


6 


0.01 ± 


0.05 


-0.06 ± 0.07 


0.0 + 0.05 


0.0 + 0.04 


0.36 


HD 140283 


5773 


3.66 


1.5 


-2.38 


8 


1 


-0.03 = 


t0.06 


-0.38 ± 0.04 


0.0 


-0.08 


0.43 


G 64-12 


6407 


4.20 


2.3 


-3.12 


4 





0.13 ± 


0.04 


-0.26 ± 0.02 
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Fig. 5. [Cr/Fe] ratios as determined from Cri and Crn lines 
under NLTE (filled symbols) and LTE (open symbols) as a 
function of elfective temperature (top) and metallicity (bottom) 
for selected stars. The NLTE abundances are computed with 
5h = 0. 



Fig. 6. Abundance ratios [Cr/Fe] as a function of metallicity. 
NLTE and LTE-based Cr abundances in metal-poor stars are 
marked with filled and open symbols. The evolutionary curves 
for [Cr/Fe] are computed with the CE model for the solar neigh- 
borhood adopting dilferent sets of SN II yields (a) - (c). 



to the LTE analysis of Cri lines. ICavrel et al.l (120041) showed 
that [Cr/Fe] ratios in giants increase from -0.5 at [Fe/H] -4 
to -0.25 at [Fe/H] x -2, and the scatter in abundance ratios 
is very small, cr = 0. 05 dex. Such a small di spers ion contra 



diets earlier results of Mc William et alj (Il995h and lRvan et al 



1996I) . who demonstrate a large spread of [Cr/Fe] abundances 



at a ny given me tallicity. The reason for this disagreement is 
that ICavrel et aL ( 2004,) investigated on ly g iant stars, whereas 



the analysis of iMc William et alJ (Il995l) and lRvan et"aD d 19961) 
are based on a compilation of their own LTE measurements in 
a mixed sample of stars and on other data from the literature. 

The mean LTE abundances in metal-poo r stars based on the 
Cm li nes are identical in our study and in Gratton & SnedenI 
(Il99lb . <[Cr/Fe]) = -0.05 + 0.04 dex and -0.04 ± 0.05 dex, 
respectively. Both of these values are consistent with the mean 



NLTE abundance derived from the Cri lines. This is not sur- 
prising, because we find minor NLTE effects in Cm for the 
range of stellar parameters investigated here and excellent 
agreement is achieved between the Cri and Cm fines under 
NLTE. 



Recently, iBonifacio et al.l (120091) performed an LTE anal- 
ysis of Cr abundances in very metal-poor dwarfs and giants, 
demonstrating the influence of multidimensional RT effects 
on Cri line formation in low-metallicity model atmospheres. 
For giants, they determine a factor of 2 smaller 3D- ID abun- 
dance corrections than for dwarfs. Although this result may ex- 
plain the systematic offset of 0.2 dex, which they find between 
evolved and turnoff' stars, it is not clear how these 3D effects 
will manifest themselves in full 3D NLTE RT calculations. 
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6. Chemical evolution model for [Cr/Fe] ratio 

6.1. Nucleosynthesis 

Calculations of stellar nucleosynthesis and Galactic chemical 
evolution indicate that the solar abundances of Fe and Cr de- 
rive for_abou^^/3_Jrom laai^^ from SNe 
II dMatteucci & Greggiol Il986l: ISamlandl Il998h . In both in- 
stances, the formation of these even -Z iron peak elements oc- 
curs in explosive burnin g of silicon ( Hix & Thielemannll 19991: 
Thielemann et al.ll2007l) . ^^Cr, which constitutes 83.8% of the 
solar Cr abundance, is formed as radioactive ^^Fe mainly in th e 
incomplete Si-burning region dUmeda & Nomotol2002li2()()5l) . 
The dominant isotope of Fe, ^^Fe (91.7% of the solar Fe abun- 
dance) is produced as radioactive ^^Ni at ve ry high tempera- 
tures, T > 4 ■ lO'^K dWooslev & Weaveilll995h . Since both ^^Cr 
and ^^Fe are formed as ff-particle nuclei, the total yields of Cr 
and Fe are almost independent of the neutron excess in the mat- 
ter undergoing explosive burning and, to a first order, of the 
metallicity of SN progenitors. 

Although the nucleosynthesis sites for Cr and Fe are iden- 
tified, theoretical stellar yields of these elements are subject 
to large uncertainties related to modelling the SN explosion. 
The amount of Fe-peak elements ejected in a SN II primar- 
ily depends on the explosion energy and on the the mass cut 
between the ejected material and the collapsed core (note that 
these parameters are related, although they a re often treated in- 
dependently, e.g. lUmeda & Nomotol (I2OO2I)). Less Cr is pro- 
duced relative to Fe for deeper mass cut ('Nakamura etaljl999l: 
Umeda & Nomoto 2002,) . The current crude estimates of the 
mass cut are based on the mass ejected in the form of ^*Ni 
and on the ratio ^^Ni/^^Ni, which are determined from SN II 
light curves and spectra. Other parameters, which affect the rel- 
ative production of Fe-group nuclei, are the degree of mixing 
and f allback to the remnant , ene rgy and geometry of explo- 
sion jNakamura et al.l Il999t lUrneda & Nomoto 2002 . 2005 ). 
Nagataki et al.l (Il997h demonstrated that an axisymmetric ex- 
plosion in SN II models enhances the production of ^^Cr 
All these parameters are interlinked and depend on preceding 
stages of stellar evolution determined by initial physical prop- 
erties of a star 

The yields of Fe-peak elements from SNe la are also rather 
uncertain. It is not clear what type of binary system evolution 
leads to an explosion, i.e. double or single degenerate scenario, 
and whether the explosion proce eds via delayed deton ation or 
fast deflagration mechanism (see Iwamoto et al.lll99"9[ and ref- 
erences therein). The production of neutron-rich isotopes in 
carbon deflagration models depends on the neutron excess 77, 
which in turn depend on the ignition density and fl ame prop- 
agation speed. For exampple. Ilwamoto et al. d 19991) found that 
the yield of a neutron-rich ^^Cr isotope in the widely-used W7 
model is very sensitive to 77 in the central region of a white 
dwarf model, where it is determined by electron captures on nu- 
clei. In this SN la model, the neutron-rich isotope of iron ^'^Fe 
is produced in the outer regions where 77 stems from CNO-cycle 
production of ''*N, i.e. the abundance of ^'^Fe is correlated with 
initial metallicity of a white dwarf. Thus, depending on SN la 
model parameters, diff'erent isotopic abundance ratios for Cr 



and Fe can be obtained. Some fine-tuning of a SN la model 
may be necessary to reproduce the solar isotopic composition. 

6.2. Model of chemical evolution 

The results we obtain for the chemical evolution of Cr for the 



solar neighborhood are based on the model of Franqois et al 
(2004), which is a revised version of the two-infall model 
presented bv IChiappini et al.l (Il997l) . This model assumes two 
episodes for the Galaxy formation: the first forms the halo on 
a short timescale of ~ 1 Gyr, the second, much longer with a 
timescale of ~ 7 Gy r, produces th e disc. The initial mass func- 
tion is adopted from lScalol ( 1986*) and the star formation law is 
the one proposed originally by Talbot & Arnett ( 1975). Further 
details on the model can be found in lFranqois et alj (|2004|) . 

The following prescriptions for Fe-peak nucleosynthesis 
are adopted. For SN la, we use the metallicity- i ndepe ndent 
yields of Cr and Fe computed by Ilwamoto et al. ( 19991) with 
the updated v ersion of Chandrasekhar mass W7 model of 
Nomoto et al] ( [l984 ). This model predicts that ~ 9% of the 



total solar Cr is in the form of ^"^Cr. Hence, we have decreased 
the yield of the neutron-rich isotope ^'^Cr by a factor of 4 in 
order to reproduce the solar abundance ratio of Cr isotopes. 
The overproduction of ^"^Cr can be avoided in other SN la 
models w ith lower centr al ignition density and neutron ex- 
cess (Iwamoto et al.lll999 ). We use three sets of massive star 



yields. The first set (a) is metallicity-independent yields from 
IWooslev & Weavej (1 1995 ) (WW95), computed with the SN II 
models with a solar composition. We increased the Cr yields 

20 Mq m odels following the rec- 



by a factor of 10 for the 10 
I 

ommendation of Francoi s et alj (120041) . who used this trick to 



avoid strong underproduction of Cr relative to Fe with their 
CE model, [Cr/Fe] ^ -1 at [Fe/H]= -4. We sugges t, how 



ever, that the Cr/Fe deficit obtained by Francoi s et al.l (120041) 



is not due to inadequate yields for low- and intermediate-mass 
SNe II. Inspecting their Tables 1 and 2, we find very low 
Cr and high Fe yields (from unknown source) for stars with 
2OM0 < M < IOOM0, which dominate Cr and Fe production 
at [Fe/H] < -4. WW95 provide yields only for M < 40Mq 
models, and their Cr/Fe production ratios for solar-metallicity 
stars with M > 15 Mq are supersolar In the second case (b), 
we adopt the metalhcity-dependent WW95 yields for Cr and 
Fe calculated for Z/Zq = 0, 10"^ 10"^ 10"', 1, and we per- 
form linear interpolation between these values. The third set 
(c ) differs from the case (b) in that for Z = we use the yields 
of iHeger & WooslevI (l2008h . computed for metal-free progeni- 
tors. 

The [Cr/Fe] versus [Fe/H] relations computed with the 
chemical evolution models (a) to (c) are compared with NLTE 
and LTE abundance ratios in Fig.|6] The theoretical [Cr/Fe] ra- 
tios are normalized to the solar abundances of Cr and Fe, pre- 
dicted by each of the models. The Cr abundance in the ISM at 
the solar system formation computed with the model (a) is 5.61 
dex, whereas the models (b) and (c) predict logscr.o - 5.68 
dex. The high abundance of Cr in the models (b) and (c) 
is compensated by a higher Fe abundance log epe.o — 7.53. 
The [Cr/Fe] trend based on LTE abundances derived from Cr i 



12 



M. Bergemann & G. Cescutti: Cr in metal-poor stars 



lines is qualitatively reproduced by the model (a) with ad- 
justed WW95 yields. [Cr/Fe] decreases with metallicity. This 
result is expected because our LTE abundances generally agree 
with the spectr oscopic abundances in metal-poor stars used by 
Franqois et a to calibrate SN yields in their model. The 



and was renorm alized to the photospheric Si abundance of 



Shi et al.l (l2008h V A few Cr i and Cr n lines give systematically 



flat [Cr/Fe] trend based on NLTE abundances of Cr is well fol- 
lowed by the model (b). Deviations of [Cr/Fe] from zero are 
due to mass and metallicity-dependence of yields in the WW95 
models. The Cr/Fe production ratios are supersolar for more 
massive stars, M > I8M0, with Z/Zq - 0.01,0.1, 1. However, 
the IMF-integrated Cr/Fe production ratios are subsolar for 
SNe II with Z = 0.01 and Z = 0.1 that is reflected in decHning 
[Cr/Fe] for metalhcities -2 < [Fe/H] < -1. At [Fe/H] > -1, 
[Cr/Fe] rises again due to supersolar Cr/Fe production from 
exploded massive stars with initial solar metallicity, Z = Z©. 
The total solar Cr abundance is attained due to dominant con- 
tribution (~ 60%) of SNe la to production of this element. We 
obtain sim ilar results using the new nucleosynthesis yields of 
iHeger & Wooslev (2008) for metal-free stars, model (c) in Fig. 
|6] The difference with model (b) can be seen only at very low 
metaflicity, [Fe/H] < -3. 

Evolutionary curves for [Cr/Fe], which are similar to the 



models (b) and (c ), were calcu l ated b y lTimmes et al. (1995); 
LSamland ('1998'); 'Ar gast et al.l (l2000l) : lOoswami & Prantzosi 
(2000); Kobavashi et al, ( I2OO6I) . However, good agreement be- 
tween our NLTE [Cr/Fe] trend and predictions of diff'erent the- 
oretical GCE models does not allow us to constrain the latter 
The reason is that different studies adopt different prescriptions 
for stellar nucleosynthesis, mixing in the ISM, recipes for in- 
flow and outflow, IMF, etc. We are not aware of any complete 
systematic study in the literature, focusing on the effect of these 
and other parameters on evolution of elemental abundan ce ra- 
tios. Studies analogous to that of iRomano et al.l (120051) . who 
investigated the effect of initial mass function and stellar life- 
times, are highly desirable. 

7. Conclusions 

The atomic level populations of Cr in the atmospheres of late- 
type stars are very sensitive to the non-local UV radiation field, 
which drives them away from the LTE distribution given by 
Saha-Boltzmann statistics. Using ID static model atmospheres, 
we find strong NLTE effects for the minority ion Cri, they 
are related to overionization from the low-excitation odd Cr i 
levels with large quantum-mechanical photoionization cross- 
sections. The NLTE abundance corrections to Cr i lines are pos- 
itive and increase with decreasing model metallicity and grav- 
ity. The number densities of excited states in Cr 11 are also mod- 
ified by non-equilibrium excitation processes, but the effect on 
abundances is negligible for dwarfs and subgiants analyzed in 
this work. 

Ionization equilibrium of Cr for the Sun and the metal-poor 
stars is satisfied, if we allow for extreme NLTE effects in Cr i by 
neglecting inelastic collisions with H i in statistical equilibrium 
calculations. For the Sun, the NLTE Cr abundance is 5.74+0.04 
dex, it is 0.05 dex lower than that computed from Cm lines. 
Both values disagree with the Cr abunda nce in C i meteorites , 
5.63 ±0.01 dex (this value was taken from lLodders et al.l(l2009b 



higher abundances. Since these lines are sensitive to microtur- 
bulence parameter, the abundance anomaly most likely reflects 
the shortcomings of our ID model atmospheres. Line forma- 
tion calculations with 3D hydrodynamical models are neces- 
sary to confirm this. 

The NLTE Cri-based abundances in metal-poor stars are 
systematically larger than that computed under LTE approach. 
The difference of 0.2 - 0.4 dex is due to substantial overioniza- 
tion of Cr I at low metallicity. The LTE abundances determined 
in this work using Cri lines are consisten t with other LTE 
studies dl ohnson 2002; Cavrel et al. 2004; Cohen etal.l 12004 



Lai et al.ll2b08l; iBonifacio et al.il2009l) . confirming that declin- 



ing [Cr/Fe] with metallicity is an artifact of LTE assumption 
in line formation calculations. The mean NLTE [Cr/Fe] ratio 
in stars with subsolar metallicity computed from Cri lines is 
([Cr/Fe]) - with the standard deviation cr - 0.06 dex. Using 
the Crn fines, we derive ([Cr/Fe]) = -0.05 + 0.04 dex. The 
finding that [Cr/Fe] remains constant down to lowest metalhc- 
ities is consistent with nucleosynthesis theory, which predicts 
that Cr and Fe are co-produced in explosive Si-burning in su- 
pernovae. 

The NLTE [Cr/Fe] trend with [Fe/H] is reproduced by 
most of the Galactic chemical evolution models, without the 
need to invoke peculiar conditions in the ISM or to ad- 
just theoretical stellar yields. Our theoretical evolution of 
[Cr/Fe] in the solar neighborhood, computed with the two- 
infall GCE model of Chiappi ni et al.1 (fl997h with SN la yields 
from Iwamotoetal. (1999) and metallicity-dependent SN II 
yields from Wooslev & Weaver (1995), is in agreement with 
the NLTE results. The model predicts that ~ 60% of the to- 
tal solar Cr and Fe are due to SNe la and the rest due to SNe 
II, the latter synthesize both elements in roughly solar propor- 
tions. The underproduction of Cr relative to Fe in SNe la is 
compensated by its overproduction in solar-metallicity SNe II. 
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